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ABSTRACT 


By the use of two separate mesh sizes instead of one in computing finite differences, an extrapolation to effective 


“‘zero”’ mesh size may be made in order to reduce the truncation error, as originally suggested by Richardson. 


This 


technique of ‘‘difference extrapolation” is applied to the estimation of individual differentials in the barotropic vor- 


ticity equation, and here corresponds to the use of ‘“‘second-order’’ finite differences. 


The trunecation-induced phase 


speed lag of the difference solution relative to the true solution is shown to be systematically reduced, especially for 


the shorter waves. 


Next, the extrapolation is applied with two separate solutions of the barotropic difference equation, 


with the result that the phase speeds are further improved, but at the expense of an amplitude distortion of about 10 


percent. 


improvement obtained, but the amplitude distortion remains for other wavelengths. 
extrapolation’’ are therefore felt to be unsuitable for routine use. 


This amplitude distortion may be removed for a particular wavelength, and a small further phase speed 


These methods of ‘‘solution 
The method of ‘‘difference extrapolation,’ how- 


ever, preserves the solution’s amplitude, and if used in conjunction with a suitable smoothing procedure should 
result in a net error reduction for those waves resolved by the mesh and retained by the smoothing. 


1. INTRODUCTION 


In the approximation of spatial derivatives by finite 
differences, a truncation error is made which in general 
depends upon the size of the finite space increment, as 
well as upon the wavelength and orientation of the con- 
tinuous field being estimated. The most widely-used 
procedure is the familiar centered space difference, which 
we may illustrate in the case of the first derivative of a 
(continuous) function f as 


of f(x+ Ax) — f(x—Ar 
? ~Af= . 7 7 ) (1) 
Or 2Ar 
where z is a typical space coordinate and Az is the mesh 
size. The error of this approximation is easily shown to 
be of the order (Ar)?, i.e., 


'T\\s research has been supported by the Geophysics Research Directorate of the Air 
Cambridge Research Center under contract No. AF 19(604)-4965., 


~ Of a ; 
e(f) =Af—s" -O( Az’), (2) 


This approximation is furthermore a consistent one in 
the sense that e(f)->0 as Ar->0. It is the purpose of this 
study to utilize such Az-dependence of the truncation 
error in order to improve the accuracy of the finite- 
difference approximation itself, using an extrapolation 
technique. The procedure to be described will be seen 
to be related to several proposed methods for the reduc- 
tion of truncation error and for smoothing. 


2. THE EXTRAPOLATION TECHNIQUE 


Since the accuracy of the difference approximation de- 
pends upon the mesh spacing, it occurs to one that from, 
say, two finite difference estimates, made with different 
grid increments, an improved estimate of the derivative 
Considering the estimation of 0f/dz for 
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might be made. 
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example, and denoting the two grid increments by Ax, 
and Arn, we have 
f(x+An,2)— f(r—4n,,+) 


(Af)i,9 : DAP ? (3) 
& 1,2 


where (Af), » denotes either of the two difference estimates 
of Of/Or. Since the error of (Af), is of the order (Az)’, 
we may extrapolate to an effective zero mesh size by form- 
ing the linear combination 
(Af), = P(4S)— (Af (4) 
; x aa 

where p=Ar,/Ar,>1. This “extrapolation to the limit” 
is illustrated in figure 1, and was evidently first suggested 
by Richardson [6]. 

In order to show that the estimate (4) is a systematic 
improvement over either (Af), or (Af), alone, we may ex- 
pand fin a Taylor series with remainder about +=0, 


Of , Arid’. Ax} O*/ 
S(at An) =f(z) + Ax; —v 2 dr? 6 Ov 
Arjo'f Ax} (O° 5 
+ 24 Ox! 30 (ae ), rey ) 
. of p’Ar? o*f 
f(t An) =f (rtp An)=f(4) +p Any +5 55 
prAat OF, pda OY ane ) (6) 
6 Ot 24 dz** 120 LOz4), 0 


where 6, . . . , 6 are suitable points in the vicinity of 
the origin. Inserting these expressions into (3) and (4), 
we find the truncation error 


!_ var’). (7) 
I 


e(N=(No-$ 

; fa) 

This result was only to be expected, however, since the 

p’-extrapolation in (4) was designed to eliminate the 
dominant (Az)? error-dependence of (Af); ». 

This extrapolation technique has been used with some 
success in steady-state problems of engineering (Salvadori 
[8]), and is a relatively well-known procedure in the 
numerical analysis of linear differential equations (see, 
for example, Buckingham [1] or Hartree [3]). In the 
nonlinear, time-dependent equations typical of dynamical 
weather prediction, however, we have no guarantee that 
its use will result in a systematic improvement. In the 
first place, the extrapolation might be performed at each 
time step in the calculation of the non-homogeneous terms, 
or might be applied to the entire solutions from the 
separate grids. Moreover, in the presence of several 
wavelengths in the field of the dependent variable, the 
extrapolation may improve the difference estimates for 
only certain waves and fail to provide an overall error 
reduction of the order implied by (7). In spite of these 


misgivings, it seems worthwhile to give the extrapolation 
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Figure 1.—The Richardson extrapolation technique for the 


reduction of truncation error, applied to the calculation of the 
derivative Of/Or by centered differences. 


technique some further consideration, particularly in view 
of its relations to other computational schemes discussed 
below. 

For the sake of clarity, let us select Az,/Az,=2, and 
again consider the estimation of 0f/Oz by centered differ- 


ences. Expanding (4) we find, with p’=4, 


(Afo=s (Api-3 (Af)o, (8) 


which, upon use of (3), in turn yields 


] 


(Afo= jan 


(f-2as—8f -ast Sfaz—foraz); (9) 
Ly 


where f_.s, denotes f(r—2Azr), ete. This expression is 
recognized as just the “five-point”? approximation to the 
first derivative. 

The extrapolation technique may also be employed in 
the estimation of other differentials, of which the familiar 
Laplacian V?f is a convenient example. Proceeding as 
before, we find 
p(A* Pi (A*D2 


a (10) 





(A?f) >= 


where (A*f),» are the difference estimates of V?/ with a 
truncation error O(Az*’). The error of the extrapolated 
estimate (A*f)y> is now O(Az*). If we select p?=2, we 
then have simply 


(A? f)o=2(A?f)1— (A? fo. (11) 
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Th.s formula is the same as that found by Knighting [4] 
in considering both the conventional five-point Laplacian 
estimate and that involving the four “corner” points at 
a distance y2 greater from the central point. This 
scheme is generated by a 45° rotation of the grid axes, 
and is an improvement on the directionally-averaged 
estimate of the Laplacian proposed by Thompson [10]. 

As a final case of meteorological interest, we note that 
the extrapolation method may be used in estimating the 
Jacobian J(a, 6) of two scalar variables a and 6. The 
conventional estimate 


](a,b)= (4.A2,7)~"[(g41, ;—@i—1, )(B;, 341 — x, 3-1) 


— (4, 541 —@s, g—-1) (D441, p—B 4-1, 9)] (12) 
may be combined with the analogous estimate with Azo, 
yielding 

I,(a,b) =[pI, (a,b) —I,(a,6)\[p?—1]=!, (13) 
since the dominant truncation error in each case is O(Az*). 
The error of the estimate J,)(a,b) is readily shown to be 
O(Ar'), and hence a systematic improvement over either 
I(a,b) or I,(a,b) alone. If we select p=2 as before, we 
find 


I,(a,b) =21, (a,b) —I,(a,b) (14) 


which has also been suggested by Knighting [4] using the 
45° “rotated” axes for J,(a,6). Knighting, Jones, and 
Hinds [5] have used both J,(a,6) and J,(a,b) as separate 
estimates of the Jacobian in numerical integration of a 
simple dynamical model, and find in general that the 
“rotated”’ axes estimate (J,(a,6)) vields a smoother field, 
and in some cases improved the prediction. From the 
present viewpoint it would have been interesting if both 
estimates were used as in (14). 


3. THE EXTRAPOLATION OF DIFFERENTIALS 
IN THE BAROTROPIC VORTICITY EQUATION 


Instead of the estimation of a single differential, the 
meteorological prediction problem involves the solution of 
a complete differential equation. The simplest of such 
equations is that of the linearized barotropic model, 


v0, 


rt 


o*'y 


’ ; : o'y 
tox 


=a (15) 


+1 +8 3 
where y is the stream function of the nondivergent flow, 
(’ the assumed constant zonal current, 8 the Rossby 
parameter, and z and ¢ the eastward spatial and time 
coordinates, respectively. As noted earlier, there are evi- 
dently two ways of applying the Richardson extrapolation 
technique: Either to improve the estimates of each dif- 
ferential of (15) separately, or to improve the solution by 
a suitable combination of two solutions found with dif- 
ferent mesh sizes. We shall consider both methods in this 
and ‘he following sections. 

The first extrapolation method, which might be termed 
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“extrapolation of differentials,” is based upon the fact 
that the truncation error of dy/dzr, d°y/dz, and d*y/dz* is 
O( Az’) when the usual centered space differences are used. 
The use of centered time differences in the solution, how- 
ever, introduces some awkwardness into the analysis by 
virtue of its inapplicability for the first time step. For- 
ward time differences, on the other hand, are unstable and 
consequently unsuitable. There remains the method of 
implicit differences (Gates [2]), and it will be used here 
in order to simplify the analysis and at the same time to 
employ a method applicable in practice. 

The essence of the implicit difference approximation is 
the use of information at both the current discrete time 
step 7 and at the next time step 7+ 1 in order to evaluate 
centered spatial differences. For the derivative dy/dr we 
thus have 


re) ] | 7 ) 
Ov mw 5 ((Afle+ (Af)r41), (16) 


where Af denotes the usual centered difference estimate 
as in (3). If both (Af), and (Af),,; are now improved by 
the extrapolation technique with p?=Ar,/Ar,=4 for 
convenience, we have 


ad ™ (24Az,) —"[8(Wn4t.24+1—Wn- 1,r+1 +Vm+1,1—Wm—1, 7) 
— (Wm+2,1+1—Wm 2,241 +Vm+2,2—Wm-2,2)], (17) 
where 
Vn, -=W(mAz,rAt), m=0,+1,..., r=0,1,.... 


In a similar manner we find the extrapolated implicit 
difference estimates 


OF (242) "[16 Yrs a.es1 + Ym t.e1— mr +Wm+1,9 
+-Vmn—1,2— 2m, 2) — (Wn +2, 241 +m —2, 941 
— Win. r+1+Vm+2,2+Wm-2,2—2bm,r)], (18) 
O°¥~(96A24)“"(32(Ym 2.041 — Winter tm t.e 
—Wm—2,74+1 t+ Vm-2,2— 2hmti.e t+ 2Vm-1,2—Vm-2,1) 
— (Wm 44,741 — 2Wmt2,24+1 + 2Wm—2, 241 —-Wm—4,741 
+ n+s,2—- 2Wm42,e+ 2Wm-2,2—Wm-s.7)], (19) 
SOY (12184) (16m s3.2414 Vm—1,7+1—2Wm, 241 
—Wmtt.e—Vm-1,2 t+ 2h, 2) — (Wmt2,041 + Vm-2, 241 
— Win, 241—Wm42,2—Wm—2,2+ Wm, r)]. (20) 


Inserting these expressions into (15) and assuming a 
solution of the form Wm,.=A(r)e"*, where a=2rdz/L 
with L the wavelength, we find after some manipulation 


(a+ ib) A(r+ 1 )—- (a—1b) A(r) = 0, 
where 
a=cos a—7, 
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b ams [¢ (cos? a+ cos a—8)+BArj ( aeeios )} (23) 


2Ar; cos a— | 


The so-called amplification matrix of (21) is given by 


a—ib 


G¢= (24) 


“7? 
a+ib 
from which we see that |@/’*=1, a characteristic of the 
(stable) implicit difference scheme. The solution of (21) 
for arbitrary t=7Aft may now be written as 

A(r)=W0G", (25) 
where V is the amplitude of the assumed initial condition 
Vn on Wve"; ie, A(O)=¥. From (25) we may then 
write the complete solution as 


Vin. = WG"? exp (i[ma+r tan~! (2ab/a?—b*)|}. (26) 


Recalling |@|)*=1 and after some further manipulation, 


we may write this solution as 


Vin. =W exp {i [ma—r tan~! (44*/4—A**)]}, (27) 


where 


as. -[ cos? a+cosa—s 
* ~ 
. SM sina{ t [ cos a—7 ] 
BL2/4r? [2 cosa—S . (28) 
‘tess *\ cosa—7 i 


a/2 


Comparing A* with the corresponding parameter 
BL2/4r? 
9 


( sin a/2 ) 
a/2 


Ai. 7 
d Ay sine U- 


occurring in the implicit difference solution of the baro- 
tropic vorticity equation without extrapolation (Gates 
{2]), we note that they differ only by the presence of the 
bracketed terms of (28). The variation of A* with L 
and Av, is shown in figure 2 for the selected values At=1 
hr., (7=20 m. see.~', and 6(45° lat.) =1.619* 107" em.7! 
sec.-' By comparison with the corresponding values of \ 
also given, we note A\*>) for all Av and L. 
The solutions with and without extrapolation may now 
be conveniently compared by writing (27) in the form 
Ym. W exp (ik | mAm—rAtCh)}, (30) 


where (is the phase speed of the numerical solution and 


k—=2rn/L is the wave number. Here C% is given by 
L _f 4* ; 
( a ar tan . a (31) 


and its variation with Az and Z is shown in figure 3, 
along with the corresponding data for Cy without extrap- 
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olation. The phase speeds are seen to be systematically 
increased for all ZL and Az, particularly in the intermediate 
cases, say 5Ar <1 <10Azr. If we compare these curves with 
the corresponding Rossby phase speed 


Cr p_bl? 

4r? 
of the analytic solution of (15), we may say that the phase 
speed error of the numerical solution is approximately 
halved by the use of the Richardson extrapolation tech- 
nique with p=2. 

The departure of the numerical phase speeds from the 
corresponding Rossby phase speeds of the continuous 
solution is due to both space and time truncation errors. 
The present extrapolation technique, we note, seeks to 
reduce only the spatially induced error, and it is therefore 
of interest to investigate its relative efficiency in compari- 


son to the limiting case of Ar—-0. From (28) we first see 


that 
: 2rAt 
lim \*= 7 c.. (33) 
Ar-0 4 


where (Cz is given by (32). Hence, from (31) we have 
lim C%= Oe ae | 2rAt Ce/L ’ (34) 
ar 2rdt 1—(wAt ('p/L)? 

and the difference between this expression and Cz may 
then be attributed to the time truncation error alone. 
In other words, the limiting speeds given by (34), and 
shown in figure 3, are the most accurate which could be 
found by reduction of the space mesh alone, and are a 
reasonable basis for comparison with the results of the 
extrapolation technique. In the case of Az,=100 km., 
for example, we see that as Z increases an increasingly 
large fraction of the total possible phase speed improve- 
ment is provided by the present extrapolation. We also 
note that for L>1500 km. (with the assumed values of 
At, LU’, and 8) the limiting phase speed need not be distin- 
guished from (',; this is to say that the time truncation 
affects mainly the shorter waves. 


4. THE EXTRAPOLATION OF SOLUTIONS OF THE 
BAROTROPIC VORTICITY EQUATION 


In the earlier discussion of the extrapolation technique, 
we noted that the method may evidently be applied either 
in the estimation of the individual differentials or in the 
improvement of the separate solutions themselves. The 
first method was considered in section 3 above, and we 
now turn our attention to the second approach, which 
may be termed “extrapolation of solutions.”’ This tech- 
nique is, in fact, that originally suggested by Richardson 
[6], who termed it the “deferred approach to the limit.” 

Continuing to use implicit differences for convenience, 
the difference equation approximating (15) is 


? This same limiting expression also results from the unextrapolated solutions give. by 
Gates [2]. 
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Ficure 2.—The variation of the parameter A* of (28) with Ar and L. Also shown by the dashed lines is the corresponding parameter (29) 


for the unextrapolated difference case. Here At=1 hr., U= 20 m.sec.—', and 8 (45°) =1.619X 10-8 em.~! sec.~! 
(AtAz?) — (Want, 41 — 2h. 41 + Wm—1,74+1— Wms Since the truncation error of both (37) and (38) is O( Az’), 
Win s—WVm—1. 2) +U (482°) —"(Vn 2 2 — Watt. + 2Vm—1.9 the appropriate extrapolation is of the form 
j 4 = 419, = ee ‘ , 
Cis Wm+2.r41 QW mn+i.r41 QYin l,r+1 Vm 2,r 1) Vin, p= p Vin . —VYn,r)(p?—1 ) a 39) 


+ B(4Ar)~ (Wns 1,2 — Vm Lert Vnti.7+17—-Vm 1,741) =0 


(35) Selecting p=As,/Ar,=2 for convenience, we then have 


The solution of this equation for the simple harmonic yi and o y (40) 
. e.e — . . . mT my T 6 m T 
initial conditions used earlier is readily shown to be 3 3 


Vn.c=W exp i | mar tan-! (OU )}}. (36) Combining the separate solutions (37) and (38) according 
ie to (40), and noting m,a;=ma, at the points of the larger 


_ mesh, we find after some manipulation 
Where \ is given by (29) and a=27Az/L as before. 


We may now write this solution specifically for our two Ving. = WA exp! ikl m Ar,—rAtCy) }, (41) 
mesh sizes Ay, and Ars as 


: 4A, ___ Where the amplitude factor A is given by 
Wmy.t Wexpd i [ mas—r tan”! (7 «| » (37) 
oo 


l ° 4r ° 4r. 
| —- —R cos 1 = 1 2 
and Aah { 17 8 co: [ tan (=F) 7 tan (oe )|} 
; _ 1 (4% Ss 
Wg, > wexpd i | meas 7 tan” (— x) )]}: (38) (42) 


and seek to combine them by the extrapolation technique. and the numerical phase speed Cy is given by 
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Ficure 3.—The variation of the numerical phase speed Ch of (31) with Ar and L, for the “extrapolation of differ- 
ences” technique. Also shown by the dashed lines is the corresponding phase speed for the unextrapolated difference case. The 
Rossby phase speed Cp of the continuous solution, and the limiting phase speed as Ar—0 are also given for comparison. Here Af=1 
hr., (7 =20 m. see.~', and 8 (45°) =1.619 10-8 em.-! see.-! 


o— L ant Av, and L is shown in table 1. We note that an amplitude 
2nrAt distortion of a few percent of the extrapolated solution 
, 4h, . 4nr, 
4 sin] 7 tan”! :) —sin} 7 tan”! : .) 
5 [ Fare 7 4—),? 
if 4 if 4r2 
4 cos] 7 tan ( .) —cos| 7 tan ( :) - ? act 
4—X,° 4— i," TABLE 1.—The variation of the r=1 amplitude distortion A of (42) 
os as a function of wavelength L and selected mesh size Ax;. Here 
(43) the second mesh Ar.=2Ar,, U=20 m. sec.', At=1 hr., and 


B(45°) = 1.619 X 10-8 cm! sec.-! 


In order to simplify the further analysis of these results, a 


let us assume that the extrapolation of the two solutions F Grid mesh An, in km. 
(40) is performed at each time step. It then suffices to (km.) j [ =F | : 

: ° . . 100 200 300 400 
examine the solution (41) for s=1, and the amplitude noe ieee ; otras i 
factor and phase speed are accordingly simplified. If 400 1.0966 |... |. FERED, eR teee a 

“ae . : 600 | 1.0189 |__- ERS 
this is not done, and the extrapolation is considered to 800 | 1.0084) 1.0267 | -2-- 
° > ° >» ,° 1000 | 1.0013 | 1.0114 | . . 
be applied after an arbitrary number of time steps 7, the 1500 | 1.0001 | 1.0016} 1.0050 | #10082 

: ° Pn . 2000 1.0000 | 1.0003 | 1.0013 1. 0027 
amplitude A will vary with the choice of r+ between the 2500 | 1.0000} 1.0001} 1.0004) — 1.0009 
ir a Ja ape de ° 3000 | = 1.0000 1.0000 1. 0001 1. 0003 
limits 1<A<5/3. With r=1, the variation of A with 4000 | 1.0000) 1.0000 | 1.0000 1, 0001 

6000 1. 0000 1.0000 | 1.0000 | 1. 0000 


' From the data presented by Gates [2], we may estimate that when 1=1000 km., caren: ae 
4n=200 km., Ar:=400 km., At=1 hr., U=20 m. sec.~!, and 8(45°) =1.619X 10") ¢m,~! *for L=1200 km 
sec.~! the maximum amplitude (5/3) would first occur for r~14. tfor ca 1600 km. 
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Ficure 4.—The variation of the numerical phase speed Cy of (43) with Ar, and L, for the “extrapolation of solutions” technique. Also 


shown by the dashed lines is the corresponding phase speed for the unextrapolated case. 


solution is also given for comparison. Here At=1 hr., U 


(41) is confined to the shorter waves, for which the larger 
truncation errors are also made. 
lengths, 4y—>Az; and A-1. 

The variation of the phase speed (43) of the extra- 
polated solution is shown in figure 4 for r=1 as above. 


For the longer wave- 


The phase speed improvement is seen to be most pro- 
nounced for the shorter waves; by comparison with 
figure 3 we note that approximately twice as much im- 
provement over the unextrapolated solution is here 
obtained as was obtained by the technique of ‘“extra- 
polation of differentials.” For Az,=100 km., waves of 
length >1000 km. are now moved with very nearly the 
correct speed. The present method of ‘“‘extrapolation of 
solutions” is therefore somewhat more effective than the 
method of section 3, at least as far as phase speeds are 
concerned. 

5S. A METHOD OF 


“OPTIMUM” EXTRAPOLATION 


While the solution extrapolation technique described 
in ‘he preceding section gives an overall improvement in 


20 m. see.~', and B(45°) 


The Rossby phase speed Cp of the continuous 
1.619 X 10-8 em.—! see.-! 


phase speed, the small amplitude distortion also introduced 
(table 1) may be an undesirable feature for some purposes 
and accordingly may limit the practical application of the 
method. The solution extrapolation of (39) is based 
upon the elimination of the dominant spatial truncation 
error O(Ar’), and is in this sense only an approxima- 
tion. It is possible to ‘‘extrapolate” the two solutions 
(37) and (38) in a slightly different way, and thereby to 
improve certain properties of the combined solution. 

Let us form a simple linear combination of the solutions 
(37) and (38), 


Yin. =m, — Wm, +5 (44) 


where @ and 6} are constants to be determined by the 
imposition of two conditions upon the solution ¥,,. As 
a first condition, let us require the amplitude of ¥;, to 
equal ¥, the amplitude of the initial conditions. This 
condition will then remove the amplitude factor A found 
in the previous method (at least for a certain wavelength, 
as discussed below). As a second condition, we may 
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, 


require that the solution ¥/7, move at the phase speed 
( N (p*l v,1— v2) (p?— 1) Ps (45) 


where Cy, and Cy» are the phase speeds of the solutions 
(37) and (38). 

formal extrapolation technique, and insures that the 
solution 7, will move at a more accurate speed than 
does either of its component solutions. As before, let 
From these condi- 


This condition is suggested by the more 


us select p=2, and note mya, 


tions we may then write (44) in the form 
vi... = WVexp { if ea : r tan (=) 
+hrtn(,™)]} 
a exp {if me 7 tan fam )|} 
—bvexp{ i[ mea 7 tan (= )}} (46) 


Setting r=1 as before 


MinQ. 


and thereby determine a and 6. 
and introducing the notation 


1 4, = 
1 
A, 3 tan x) (47) 
] 4dr. 
P 1 = 
A, 3 ban (=) (48) 


for convenience, we find upon equating real and imaginary 
parts of each side of (46), 


sin [4(4, —6,)| 


; ’ (49) 
sin |3(0,;— 6.) | 


sin [6,—62] ; (50) 
sin [3(@,—#@,)| 

For a selected Ar; and Ar.(=2Az,), the values of \; and 
\» vary with the wavelength L (see fig. 1). Hence 6, and 
#, are likewise wavelength dependent, and it would appear 
that 1 must be specified for a unique determination of a 
and 6. Fortunately, however, the values of a and 6 are 
not very sensitive to the value of L selected, as shown in 
table 2 below for the case Ar, =100 km. For the longer 
wavelengths we notice that a . and b -? corresponding 


to the extrapolation technique considered earlier in sec- 


The variation of the weighting coefficients a and b of (44) 


100 km., U=20 m. sec.—'!, p=2 
1 


TABLE 2. 
with wavelength, for the case Ax, 
At=1 hr., B (45°) 1.619 10-8 cm. sec. 


’ 


Wavelength ZL (km.) 


] l 
400 600 800 | 1000 1500 2000 3000 6000 
a | 1.252) 1.318] 1.330] 1.332] 1.333] 1.333] 1.333 1. 333 
b 0.358 | 0.338 | 0.334] 0.333 | 0.333 | 0.333 | 0.3383 | 0.338 
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tion 4. Since it is the shorter waves’ phase speeds which 
require the most improvement, we may select L=400 
km., the smallest permissible value (the resolution of the 
Ar, mesh). In this ease, the phase speed of the solution 


Wnts 
eo > 4r. , 
( == — -1 Th wt e. : 
w2rat L3 (= eae GAs)} (51) 


varies as shown in figure 5. Here the corresponding 
phase speeds for the methods of sections 3 and 4, as well 
as that for the unextrapolated solution, are shown for 
comparison. We note that the use of @=1.252 and 
b=0.358 in (44) has given a small improvement over the 
use of a=4/3, and 6=1/3 as in (40), which in turn is 
somewhat superior to the “extrapolation of differentials” 
of (27). This selection of a and 6 has, in effect, “tuned” 
the weighting scheme to give the greatest improvement 
at 1=400 km., and is in this sense an ‘‘optimum ’ scheme. 
For other Ar,, and a corresponding selection of @ and } 
from table 2, generally similar results are obtained. 

With this “tuning” to L=400 km., the amplitude is 
unity for this wavelength but is generally not unity for 
other wavelengths, as shown in table 3. This amplitude 
distortion, especially for the longer waves, greatly reduces 
the attractiveness of this “optimum” method, and is far 
more serious than that of the simpler method of ‘“extra- 
polation of solutions’ (39). The small phase speed im- 
provement given by this method over that of (39) would 
therefore not be a useful improvement in practical appli- 
cation. 


6. THE EXTRAPOLATION TECHNIQUE AND 
SMOOTHING 


The basic extrapolation scheme (39) attempts to im- 
prove the accuracy of the finite-difference solutions by 
the removal of some of the spatial truncation error. 
Smoothing procedures (Shuman [9]), on the other hand, 
are designed to suppress the shorter waves (for which the 
larger truncation error is made) by the deliberate intro- 
duction of additional truncation, and in this sense the 
two procedures are opposed. If we write (39) for a 
discrete variable g,, in the form 


Pm, — Pm, =(1—p*) om, (52) 


where the notations are as before, and write Shuman’s 
“three-point”? smoothing operator for the same variable 


The variation of the amplitude with wavelength of the 
400 km. in order to 


TABLE 3. 
extrapolated solution (46), for the choice L 





determine a (=1.252) and b (=0.358) in (49), (50). Here Ax,= 
100 km., and p, U, At, and B are as in table 2. 
Wavelength LZ (km.) 
tis | eee | a ie 
400 600 200 1000 | 1500 2000 3000 4000 000 
1.000) 0.915 | 0.899] 0.895 | 0.894 0. 894 0. 894 0.894 | 0 804 
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FicurE 5.—The variation of the numerical phase speed Cy of (51) with L for the “optimum” solution extrapolation at L=400 km., for 


the case Ar;=100 km. 


- ? 77 . P . 
Also shown for comparison are the corresponding phase speeds Cy (43) for the usual solution extrapolation 


technique. CX (31) for the difference extrapolation technique, Cy for the unextrapolated solution, and Cp (32), the Rossby (contin- 


uous) phase speed. Here At=1 hr., U 


in the form 


i i 
i+ i Em+i1 7 em-1 ws 
Ym 1— woh, —_ (1 —— 955 ( > -1), (53) 


we may examine this relationship in more detail. 

In (53), the superscript 7 denotes an ‘‘unsmoothed” 
value, and 7+ 1 denotes a (once) “smoothed” value, while 
smoothing index”. We notice that on a formal 
basis the larger-mesh solution ¢gm2 corresponds to the local 
smoothed variable of Shuman, the smaller-mesh solution 
¢m corresponds to the local unsmoothed value, and the 
extrapolated solution ¢/, corresponds to the local space 
average of the unsmoothed variable of Shuman. With 
this interpretation the two procedures complement each 
other, and might well be employed together. For ex- 
ample, one could apply an extrapolation technique and 
then apply a smoothing operation, with the result that 
the behavior of those wavelengths retained by the smooth- 
ine would be better than if smoothing alone were applied. 
With Av,=100 km., say, and p=2, this could be accom- 
plished with (40) and a smoothing operator of the Shuman 
type designed to suppress waves of length <600 km. 


5S4984— 61 2 


ae 


m is the 


20 m. sec.~!, and 8(45°) 


1.619 * 10-8 em.~! sec.~! 


7. CONCLUSIONS 


The technique of extrapolating finite-difference esti- 
mates to effective zero mesh size systematically reduces 
the truncation error of solutions of the linear barotropic 
vorticity difference equation. This technique increases 
the numerical phase speeds toward the continuous solu- 
tion’s speeds, especially for the shorter waves, while 
preserving the solution’s amplitude. The extrapolation 
of entire difference solutions found with two different 
mesh sizes results in an even greater phase speed improve- 
ment, and may be adjusted to give optimum results in the 
vicinity of selected wavelengths. While this ‘‘extrapola- 
tion of solution’ technique improves the solution’s phase 
speed, it introduces an amplitude distortion of certain 
wavelengths. The ordinary solution extrapolation method 
introduces about a 10-percent amplitude increase for the 
shorter waves only, while the ‘‘optimum” method intro- 
duces about a 10-percent amplitude decrease for all of the 


“ 


longer waves, and is consequently unsuitable for routine 
8 | s 

use. The method of “extrapolation by differentials,” for- 
mally equivalent to the use of the higher-order difference 





a 


approximations, may then be the best that can be expected 
from the extrapolation Whether the in- 
creased computation required by even this method is jus- 
tified by the increased accuracy must await an actual 
numerical integration. Such a test is particularly impor- 
tant in view of the general observation that higher-order 
difference schemes often fail to provide the accuracy 
expected (Richtmyer [7}). 

Although the implicit difference scheme has been em- 
ployed in this analysis for convenience, the techniques 
described are in no way restricted to it. Noting that the 


techniques. 


numerical phase speeds of solutions found with the implicit 


difference scheme are less accurate than those found, say, 
with the first-forward-then-centered scheme (Gates [2}), 
an extrapolation technique with the latter difference 
scheme could be expected to give correspondingly more 
accurate results. 
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ABSTRACT 


Winds and temperatures to above 80 km. measured during the International Geophysical Year by grenade, 
sphere accelerometer, and Pitot-static tube experiments at Churchill, Manitoba are combined with radiosonde data 
The resulting analyses are discussed in reference to 10-mb. constant pressure charts showing the 
Periods in July, August, and December 1957 and January-February 


in time sections. 
horizontal circulation near the 30-km. level. 
1958 are covered. 

The nature of the circulation in successive layers of the stratosphere and mesosphere is suggested by this inves- 
tigation. In summer, irregular wind flow that is predominantly from the east appears in the middle and upper 
stratosphere. Strong easterlies in the lower and middle mesosphere are bounded in the upper mesosphere above 
75 km. by a sharp vertical wind shear layer with highly variable west winds above, suggesting the existence of rapidly 
moving, intense cellular circulations near the mesopause. In winter when the boundary of the polar night is close 
to Churchill, such circulations seem to extend into lower layers until in January they occupy the entire mesosphere. 


Strong west winds develop in the mid-stratosphere of wintertime, and a temperature minimum appears near 


the 30-km. level. 
by sudden warming of the entire stratosphere. 


During January 1958 the strong westerlies and the temperature minimum as well were disrupted 
The vertical structure of this warming effect, as it moved westward 


across the Atlantic and first appeared at the 40-km. level over Churchill, is shown by means of a cross section. 


1. INTRODUCTION 


Until recent vears there has been a tendency to consider 
conditions from 20 to 90 km. above the earth either as 
quiescent or else dominated by zonal motions that change 


little from day to day but substantially from season to 


There has been great uncertainty about the 
magnitude of interdiurnal variations at these high levels. 
Data of the sort presented here help to shed light on this 
problem. 


season. 


A remarkable series of temperature and wind observa- 
tions was made at Churchill, Manitoba during the Inter- 
national Geophysical Year by means of grenades, sphere 
Pitot-static tubes aloft in 
The grenade experiments are described, 
data presented, and copious references given by Stroud, 
Nordberg, and Walsh [15] and Stroud, Nordberg, Bandeen, 
Bartman, and Titus [13, 14], while accelerometer and 
static tube experiments are similarly described by Jones, 
Peterson, Schaefer, and Schulte [5, 6] and Ainsworth, Fox, 
and LaGow [1], respectively. 

In this paper, the results of these experiments are pre- 
sented in a somewhat different manner to bring out 
further details of the temperature distribution and cir- 
culation of the upper stratosphere and mesosphere. In 


accelerometers, and carried 


rocket vehicles. 


each: of the four cases discussed, radiosonde data are 
used to determine the analysis up to the 30-km. level 
during a two- or three-week period centered on the days 


of the rocket soundings. For higher altitudes, data from 
rocket experiments are available for two or more obser- 
Many facets of 
the data and their meteorological implications already 


vations during each of the four periods. 


have been discussed in papers by those responsible for 
experiments. | Unavoidably, 
repeated here in developing additional conclusions. 


these some points are 

The reliability of conclusions drawn in the following 
discussion is determined by the accuracy of the data. 
Stroud et al. [14] state that the errors of the grenade 
experiment data below 75 km. are generally less than 
+3 C., ' +15 
degrees for wind speed and direction, respectively. Above 


for temperature, and +5 m.see.~“' and 
75 km. the size of the errors may increase rapidly by a 
factor of as much as ten at the very highest points. Jones 
et al. [5, 6] estimate the probable temperature error of the 
accelerometer experiment to be about 2 percent (4°—-6° C.) 
below 75 km. and 5 percent (10°-15° C.) above 75 km. 
Furthermore, the dependence of the acceleration upon 
air density and drag coefficient limits the useful range of 
the experiment to the layer from 25 km. to 80 km. 
Temperatures from the static tube experiment are derived 
from density data having an estimated measurement error 
of less than 2 percent at 80 km. (Ainsworth et al. [1]). 
The wind and temperature values given for the rocket- 
sonde observations are averages for the layer between 
points at which two grenades exploded; thus, maximum 
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JULY 1957 


Figure 1. 


and 0530 amv July 24, and a Pitot-static tube experiment at 2200 cmv July 29. 
Grenade temperatures and winds and selected winds from special radiosonde 


data below 30 km. and on grenade data at higher levels. 
observations are included. 
indicate cold and warm areas of diagram, respectively. 


Time section for Churchill, Manitoba, July 17-29, 1957, including data from rocket grenade experiments at 0416 cmt July 2 


Fronts given by dotted lines, tropopause by heavy line, thermal boundary by dash-dot line. 
Barbs representing wind speed are placed at the tail end of the wind shafts 


22 


Temperature analysis in °C. is based on radiosonde 


C and W 


which point in direction toward which wind is blowing, with north at top of diagram and east at the right. 


and minimum values tend to be smoothed out to some 
extent. Furthermore, random errors will have opposite 
signs in adjacent layers when the source of error is diffi- 
culty in accurate timing of sound arrivals from the 
grenade at the boundary between the layers. In exces- 
sively thin layers, less than 3 km., errors may be greatly 


magnified. Where highly incompatible values suggested 


the effect minimized by combining 


adjacent lavers. 


such errors, was 


2. THE CASES OF JULY AND AUGUST 1957 


In their broad features, the time sections for July (fig. 1) 
and August (fig. 2) 1957 are remarkably similar. In the 








51 


he 














Aprit 1961 MONTHLY WEATHER REVIEW 127 
| T T T | T i T | | T 1 ] ] T | | ] qT 

90+ %, a 

92 
r * ROCKET GRENADE DATA SM 1.04 SM 1.05 q out 2.08 

@ RADIOSONDE DATA “ee z 
80}- “w somrs hte Cc M 

“TV 25 mPs 4 


























Km 
































o—. 
> 















es eee” 
— 
——— = \ 






























































16 


Figure 2.—Time section for Churchill, August 8-28, 1957. 


lower stratosphere, in each case, there was marked counter- 
action to tropospheric activity; e.g., warming occurred in 
this layer simultaneously with the tropospheric cooling 
following cold front passages. Vertical wind shear, operat- 
ing in the troposphere to produce a wind maximum at the 
tropopause, was reversed within this layer so that light 
winds were found at the top of the laver. There, near 18 
kn., we find a poorly defined surface of minimum tempera- 
ture having elevation changes out of phase and tempera- 
ture changes in phase with those of the tropopause. This 
boundary seemed to mark the upper limit of tropospheric 
inf'uence upon stratospheric conditions, for in the 20- to 


17 
AUGUST 1957 





18 


Rocket grenade experiments are for 1600 cmt August 12, 0230 cmt August 
20, and 1408 cmt August 25. 


Explanation as in figure 1. 


30-km. layer the winds remained light and temperature 
constant from day to day. During the July period and 
continuing into the August period, the temperature varied 
by only a few degrees, suggesting that there was little 
vertical motion and almost perfect radiative equilibrium 
in the laver. 

From 18 km. up to the stratopause, the thermal maxi- 
mum at the top of the stratosphere near 50 km., the 
temperature increased from about — 50° C. to about 0° C, 
The double temperature maximum at the stratopause on 
July 29 and again on August 12, already noted by Ains- 
worth et al. [1], is an interesting phenomenon that might 
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10-mb. chart for 1200 Gm, July 25, 1957 (from [18]). 
160-m. intervals; intermediate contours given by 
*. inter- 


Figure 3 
Contours at 
long-dashed lines; isotherms given by dashed lines at 6° ( 


vals. H is for high pressure center, L for low pressure center. 

















10-mb chart for 1200 Gut, August 15, 1957 (from [18]). 
explanation as in figure 3. 


FIGuRE 4. 


be due either to differential heating of layers or to advee- 
tive processes. In the mesosphere above 50 km., the 
temperature decreased slightly up to 60 km. but then 
decreased at an average rate of 4° C. km.~' up to about 


Figure 5.—10-mb. chart for 1200 amt, August 25, 1957 (from [18)). 
Explantation as in figure 3. Colder air is in the Tropies and near 
low center in Canada. 


80 km. Although on July 24 and again on August 25 a 
surface of minimum temperature defining the mesopause 
or top boundary of the mesosphere appeared near 80 km., 
the rate of temperature rise above that level was very 
slight. 

During July, winds at 10 mb. (30 km.) over North 
America were everywhere from the east as illustrated by 
contours of figure 3. On July 22 and 24 at Churchill 
(fig. 1), the east wind tended to increase with height above 
30 km. to speeds exceeding 50 m. sec.~! in the middle and 
upper mesosphere where the wind backed to northeast or 
northerly. If the temperature variation between the two 
observations and the isotherm gradients indicated by the 
wind shears are taken at face value and interpreted with 
the aid of the thermal wind equation, the conclusions 
made in the following paragraphs are valid. 

The air in the 40- to 50-km. layer cooled by about 7° C. 
in the 2-day period. The change of vertical wind shear 
from easterlies diminishing with height to easterlies in- 
creasing with height requires that greater cooling took 
place to the south of Churchill than to the north. 

In the layer from about 45 to 64 km., the observed 
vertical wind shear of about 2 m. see.~' km.~! would have 
required a north-to-south temperature gradient of 7° C. 
in 10° of latitude, although this rate of horizontal tem- 
perature variation would have had to extend beyond the 
vicinity of Churchill only if vertical wind shear did also. 
However, since the actual wind was nearly parallel to 
the thermal wind, the advective wind was negligible and 
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litle horizontal temperature advection or vertical motion 
along isentropic surfaces is indicated. 

At 83 km. on July 22, a remarkably low layer-mean 
temperature of —112° C. was reported; the decrease of 
easterly wind speed reported between 80 and 83 km. 
would have permitted even lower temperatures to the 
north of Churchill. 

By the end of August (fig. 2) temperatures in the 
troposphere and middle stratosphere were about 5° C. 
cooler than in July. Cooling was more marked at the 
stratopause where the temperature was more than 10° C. 
lower on August 25 than on July 22. The north-to-south 
temperature gradient that was present in July in the 
upper stratosphere near Churchill seems to have been 
destroved by August 25 when the winds in that layer 
showed no definite vertical shear. This apparently 
signals the beginning of that part of the year when 
temperatures near the level of the stratopause decrease 
toward the north. Indeed, the 10-mb. charts for August 
15 and 25 (figs. 4 and 5) show that such seasonal tem- 
perature and wind changes were taking place in the middle 
stratosphere with gradual formation of a cold Low over 
northwestern Canada. The 10-mb. temperatures for the 
remainder of the year at Churchill continued to decrease, 
particularly rapidly in September and October. Westerly 
winds appeared at 10 mb. early in September and increased 
in strength throughout the autumn. 

In the mesosphere up to 65 km. during August, a 
persisting north-to-south temperature gradient is indicated 
by the increase of the easterly wind component with 
height. Although the August 12 sounding was cut off 
near 65 km., it shows easterly winds as strong as those 
reported in July. However, in the 50-75-km. layer of 
the mesosphere, the two later August rocket soundings 
show progressively lighter easterly components and lower 
temperatures. 

Near the mesopause there was a pronounced vertical 
wind shear with wind shifting from strong easterly to 
strong westerly, indicating a strong horizontal temperature 
gradient of about 2° C. per latitude degree toward the 
north in the 70-S0-km. layer on August 20 and 4° C. per 
latitude degree in the 78-83-km. layer on August 25. 
The temperature of the mesopause was progressively 
warmer in each observation from July 22 to August 25. 
However, the horizontal temperature gradient required by 
the vertical wind shear across the mesopause on August 
25 would have permitted the —112° C. temperature of 
the July case to be found only 5° latitude to the north of 
Churchill. 

Murgatroyd [9] presented a vertical temperature dis- 
tribution for summer at 60° N. that is confirmed in most 
respects by the Churchill summer soundings. Among the 
principal points of difference are his relatively high 
temperatures of +21° C. and —93° C. at the stratopause 
and mesopause respectively. His temperatures for the 
mesosphere between these levels do not differ greatly 
fron: those of the grenade experiments since he gives a 
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greater temperature lapse rate immediately above the 
stratopause and a lesser lapse rate just below the meso- 
pause than do the grenade experiments. The lapse rates 
of the grenade experiments thus favor a maximum of 
vertical mixing in the layers just below the mesopause. 
The winds observed differ in two respects from those 
given by Murgatroyd for summer conditions near 60° N. 
The Churchill data indicate higher wind speeds, exceeding 
25 m. sec.~' in the 60- to 75-km. layer, and a shift to west 
winds near SO km. elevation instead of near 90 km. 


3. THE CASE OF DECEMBER 1957 

The time section for December 6-20, 1957 (fig. 6) 
shows a troposphere 25° C. colder than it was in summer. 
Above the tropopause, with average height about 1 km. 
less than in August, the lower stratosphere had cooled 
only 10° C. The middle stratosphere underwent more 
remarkable changes, the wind having increased to westerly 
25 to 50 m. sec.~' as a 30° C. decrease of temperature at 
28 km. produced an inversion base at that altitude. A 
broad trough with coldest air at the pole dominated 
Canada at the 30-km. level (fig. 7). 

Like the lower stratosphere the stratopause cooled only 
10° C. The position of the stratopause at 51 km. is at 
variance with the position near 60 km. given by Murga- 
troyd [9] for winter at this latitude. In fact, these experi- 
ments show a remarkable temperature minimum near 60 
km. The height of the stratopause in January (fig. 8) 
was indefinite but closer to 60 km. The 
extreme temperature changes above 60 km. during the 
December 12-14 period are in strong contrast to the 
regularity in temperature patterns of the summer meso- 
Coupled with winds in the neighborhood of 100 


somewhat 


sphere. 
m. sec.~', these temperature changes suggest the existence 
of violent storms in the winter mesosphere, albeit they are 
tempests in a near vacuum. Determination of the geo- 
graphical extent and speed of motion of these storms 
requires a synoptic network of rocketsonde stations such 
as that now being activated over North America. 

The wintertime temperatures of the upper mesosphere 
at Churchill were much warmer than those measured 
there in summer or at lower latitudes in winter. Although 
warming in the stratosphere and mesosphere is generally 
explainable by a reasonable amount of sinking and com- 
pression of air from higher layers, Kellogg {7} has described 
a process by which heat may be liberated at or above the 
through recombination of 
brought down from ionospheric layers rich in this element. 


mesopause atomic oxygen 


4. THE CASE OF JANUARY-FEBRUARY 1958 


Even before the rocketsonde data became generally 
available, the radical changes in stratospheric circulation 
and accompanying explosive warming of the stratosphere 
during the weeks following mid-January 1958 had been 
the subject of several papers (Boville [3]; Palmer [10]; 
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DECEMBER 1957 


Time section for Churchill, December 6-20, 1957. Rocket grenade experiments were at 0400 Gmr December 12 and 2100 
am? December 14. Explanation as in figure 1. 


Figure 6 














Scherhag [{11, 12]; Teweles and Finger [16]; Teweles, 
Rothenberg, and Finger {17]). 

Because most frontal activity in the troposphere was 
displaced far to the south of Churchill during the period 
January 22 to February 6, only relatively minor changes 
in tropopause height and temperature are revealed in 
figure 8. 


In the lower stratosphere, unusual warming 


after January 26 culminated in a temperature of —41° C. 
at 100 mb. on February 3. In mid-stratosphere by 
January 25, the 30-km. temperature had dropped below 
—70° C., and by January 27 the altitude of the temper- 
ature minimum was raised to 33 km. This temperature 
minimum was depressed to 25 km. in the succeeding days 
as strong warming occurred in the upper stratosphere. 
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10-mb chart for 1200 Gmt, December 15, 1957 (from 
{18]). Explanation as in figure 3. 


FIGURE 7. 


For the period January 25 to 29, an unprecedented con- 
centration of rocket sounding data is available for recon- 
structing conditions at high levels. Accelerometer spheres 
launched by Nike-Cajun rockets on January 25 and 29 
gave both up-leg and down-leg information. Down-leg 
information was given by a similar sphere, borne by one 
of two grenade-carrying Aerobee rockets launched on 
January 27. The temperatures derived from 
accelerometer experiments become consistent with radio- 


these 


sonde data at and above the 25-km. level, and above the 
ceiling of the radiosonde data remain consistent internally 
and with grenade data up to 65 km. 

Warming took place slowly at Churchill between 40 and 
45 km. from January 25 to 27, then rapidly between 30 
and 60 km. from January 27 to 29. The 4-day temper- 
ature increase of almost 70° C. in the 38- to 41-km. layer 
at Churchill is greater than any ever recorded at lower 
levels. 

The circulation and thermal structure at levels above 
30 km. on January 27 can be reconstructed by hodograph 
analysis (fig. 9) of grenade experiment winds for that day 
(shown on fig. 8). From the tropopause up to 65 km., 
there is remarkable agreement in the winds from these 
two experiments. A southwest wind of 30 m. sec.~! was 
observed just below the tropopause at 8 km. as a trough 
approached from the west. The wind decreased upward 
to a minimum near 14 km. and increased again to a 
maximum of about 75 m. sec.~' from the north-northeast 
at 33 km. The resulting shear vector in the 8- to 33-km, 
layer placed the warmest air toward the northwest. Thus 
there was a cold Low to the east and a warm High to the 
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west of Churchill with maximum intensity at 33 km. 
There have been few cases in which the level of a polar 
night jet stream of such strength has been so clearly 
defined. Above the 33-km. level, the decrease in strength 
of the jet required a reversal of temperature gradient. 
Thus the Low became a warm Low and the High became 
a cold High, both systems decreasing in strength with 
elevation. These features characterize the 33-km. level 
as a null level, as described by Faust and Attmannspacher 
[4] (or see Lamb [8]}). 

Near the 45-km. level the wind shifted around through 
easterly to become strong southerly at 52 km. Apparently 
the vertical axis of low pressure sloped from east to west 
across the Churchill area at about the 45-km. level, with 
the low center itself south of Churchill at that level. The 
south wind increasing with height in the 45- to 52-km. 
laver required a warm High to the east and a cold Low 
to the west, both systems increasing in strength with 
elevation to at least 52 km. The edvective wind in the 
33- to 52-km. laver can account for a warming of about 
30° C. in two days. Since a warming of 60° C. was 
actually observed at 40 km. during the two days ending 
on January 29, and since the final temperatures are 
undoubtedly anomalously high for that level, subsidence 
as well as advection must be called upon to explain such 
a great warming. 

Above 52 km. the circulation apparently retained its 
general form and strength to at least 75 km. Although 
the reported winds were very erratic, mean winds for 
overlapping 10-km. layers up to 75 km. were quite 
consistently from the south. 

As in the December case, the alternation of warm and 
cold air layers was in sharp contrast to the steep lapse 
rate of the summer mesosphere, and the temperature in 
the upper portion of the mesosphere was much warmer 
than in summer. 

In the 0604 Gur sounding, a remarkably strong west 
wind was reported along with high temperature at 79 km. 
By 1845 emt, west winds were reported at all levels above 
77 km. Expansion of the warm layer near 80 km. during 
the next two days is indicated by the falling-sphere data. 
Here we see the approach to Churchill of a system similar 
to that which moved away from Churchill between 
December 12 and 14. The reality of the large variations 
that are reported in wind speed and direction from one 
level to the next or between observations only a few hours 
apart must be more firmly established before great 
confidence can be placed in the features of the circulation 
inferred from the observations. Even so, it is now certain 
that the wintertime circulation of the mesosphere, in the 
subarctic, is very active and thus of very special interest 
to the meteorologist. 


5. EVENTS IN SURROUNDING AREAS 
AT TIME OF WARMING AT CHURCHILL 


The 1200 emt 10-mb. charts for January 25 (fig. 10), U.S. 
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Figure 8.——-Time-section for Churchill, January 22-February 6, 1958. Rocket grenade experiments were at 0604 GmT and 1845 GMT 
January 27. Sphere accelerometer experiments were at 1912 Gmr January 25, 1845 cmv January 27, and 1906 GmT January 29. 


= 


Explanation as in figure 1. 


Weather Bureau [18], and for January 29 (fig. 11), Behr 
et al. [2], help to place the Churchill rocket sounding data 
near the 30-km. level in perspective with events that were 
taking place in surrounding portions of the hemisphere. 
These charts confirm the existence of the relatively cold 
Low revealed by the wind and its vertical variation near 
this level on January 27. In the 4-day interval, this Low 
rapidly approached Churchill from a position over southern 


Greenland. Meanwhile, the warm High in the Gulf of 


Alaska moved poleward across Alaska into the Aretic 
Ocean. 

During early January, a section of the polar night jet 
stream, which was nearly circumpolar in December (fig. 7), 
slowly arched poleward across Alaska. By January 25, 
the axis of this current stretched across the North Pole 
southward toward Churchill, and the Arctic was being 
ventilated with warm air from sunlit latitudes. The 
exceedingly cold air that formerly occupied the Arctic 
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FIGURE 9. 
tail at origin; arrow points in direction wind is blowing. 
height of mid-point of layer to the nearest km. 


Two or three points were averaged where indications were indefinite. 


Ilodographs of the wind observations plotted in figure 8. 
Circled dots represent grenade experiment winds, numbers at dots give 
Arrows in diagram give vertical wind shear in layers having relatively uniform shear. 


Radiosonde winds given by dots to represent head of arrow with 


Shear arrows for layers are superimposed below hodograph. 


Relatively warm direction in the shear layer is indicated by W, cold direction by C; thus, W/C indicates the direction of greatest 


static stability with relatively warm layer lying above relatively cold layer. 
Shear of the highest layers, given by dashed arrow, is somewhat indeterminate. 


shear lying between layers with shear. 


split into two cells, one associated with low pressure over 
Siberia (see Behr et al. [2]), the other centered in the 
western portion of the low pressure area approaching 
Churchill from the east. 

Hodograph analysis for January 27 (fig. 9) when applied 
to the 10-mb. circulation pattern leads to several inter- 
esting conclusions. The polar night jet stream at Churchill 
had its maximum strength at 33 km. above which it lost 
strength. Further interpretation is suggested by the pattern 
of contours and isotherms found at 10 mb. on January 25. 
The axis of the 10-mb. Low center, located just southeast 
of Hludson Bay on January 25, from hydrostatic consider- 
ations must have sloped with elevation in the direction 


of coldest air. If the axis of coldest air also sloped in the 


Near 60 km. there is a neutral layer without significant 


same direction and the contour-isotherm relationship did 
not change with height, the vertical wind structure at 
Churchill can be explained quite convincingly. At 33, kin. 
where the wind shear direction reversed, the cold air was 
centered directly over Churchill, but the low pressure cen- 
ter was still to the east. Above 33 km., the cold air was 
centered west of Churchill to at least 51 km., but northerly 
winds show that the low center remained to the east and 
southeast up to 45 km. At this level the Low in its 
westward movement passed south of Churchill between 
observations. The uniformly strong wind shear between 
37 and 52 km. showed the presence of a tight temperature 
gradient with warmest air toward the east, like that in 
the southeast quadrant of the Low at 10 mb. (figs. 10 
and 11). 
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Ficure 10.—10-mb. chart for 1200 Gut, January 25, 1958 (from 
{18]). Explanation as in figure 3. 














10-mb. chart for 1200 Gat, January 29, 1958 (data 
from [2]). Explanation as in figure 3. 


Figure 11. 


The 4-day temperature change at 10 mb., computed 
from the isotherm patterns of January 25 and 29 (fig. 
12), shows a remarkable rise of more than 54° C. in the 
region Greenland and Labrador. The even 
larger temperature change at the 40-km. level above 
Churchill during the same period was most certainly 


between 


Figure 12.—Temperature change in °C. at 10mb. in the 4-day 
period January 25-29, 1958, obtained by graphical subtraction 
of isotherm patterns in figures 10 and 11. W indicates area of 
warming, C area of cooling. 


connected with the warming farther east at 10 mb. This 
gives convincing proof of the tremendous lateral shift with 
altitude of stratospheric patterns. In this case, the 
leading edge of the warming sloped upward 10 km. in a 
horizontal distance of 2,000 km., the same 1:200 slope 
typical of warm front surfaces in the lower troposphere. 
If the wind reversal near 45 km. is taken as evidence 
that the low-center position was very close to Churchill at 
that level, and if the low-center position at 30 km. is taken 
from the January 27 10-mb. charts of Behr et al. [2], then 
the axis of the low center tilted westward 1,500 km. in a 
vertical layer 15 km. deep, a slope of 1:100. The 1:200 
slope of the warm air indicates that it was overspreading 
the top of the low center to diminish the intensity of the 
center first near 40 km. and gradually at lower levels. 
The net effect is demonstrated by the great loss of kinetic 
energy at 30 km. over North America during the follow- 


ing week and over the remainder of the hemisphere by mid- 
February (see 10-mb. charts in Behr et al. [2] or U.S. 


Weather Bureau [18]). Simultaneously, the general level 
of the 10-mb. surface in the Arctic was substantially 
raised, 

A cross-section from Naknek, Alaska through Churchill 
to Berlin, Germany on January 29 (fig. 13) illustrates the 
great horizontal and vertical extent of the warming at the 
time of the 40-km. warming at Churchill. Cooling with 
the advance of the cold sliver near 30 km. in the lower for- 
ward portion of the system is almost as impressive as the 
warming in the upper rear portion. While the cooling 
did not produce temperatures that were anomalous for the 
entire hemisphere as did the concomitant warming, it cid 
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Fictre 13.—Cross section, January 29, 1958, near 55° N. from Naknek, Alaska (station 326) through Churchill (station 913) to Berlin 


Germany (station 384). Explanation generally as in figure 1. 


entirely on sphere accelerometer experiment AM 6.03 for 1906 amr. 


Temperature analysis at high levels in vicinity of Churchill is based 


Identification numbers of radiosonde stations from which data 


were used are shown along bottom of chart at projected location of station. 


bring anomalously low temperatures into the middle 
latitudes. All indications are that the stratospheric 
features in the cross-section were slowly moving toward 
the west. Thus, the warmest air at 25 km. crossed over 
Berlin on January 25, passed to the south of Greenland 
on January 29, as shown in figure 12, and then diffused 
across North America in the first three days of February. 

Figure 8 shows the upwelling of the coldest air to above 
30 km. that immediately preceded the intensive warming 
near the 40-km. level. The core of the polar night jet 
stream seems to have been wedged between these layers 
of climactic cooling below and incipient warming above. 

The interaction of vertical motions and jet stream 
accelerations around the vertically tilting core of low 
pressure must account for the observed temperature and 
pressure changes. Two-dimensional sections do not reveal 
the exact linkage by which these motions are maintained 
for long periods in a tremendous volume of atmosphere. 
However, increasingly detailed descriptions of these phe- 
hoiena may eventually lead to a successful mathematical 
formulation of the problem. 


6. CONCLUSION 
The temperature and wind observations utilized in this 
paper were among the first to become available to great 
heizhts at such high latitudes. The large amount of 


information deduced from this limited number of sound- 
ings is a strong recommendation for gathering much more 
temperature and wind information to the highest possible 
levels. 

That there are several lavers of the atmosphere having 
their own characteristic types of horizontal circulation is 
brought out by the charts displayed in this paper. Un- 
doubtedly, the dominating cause is the division of the 
atmosphere into layers of radically different static sta- 
bility by the variable effectiveness of its constituents in 
absorbing radiation at different altitudes, latitudes, and 
seasons. However, the circulation at any level is tied 
hydrostatically to the circulation at adjacent levels by the 
temperature distribution in the intervening layers. Feed- 
back of energy from one level to another may be great or 
small depending upon the static stability, baroclinicity, 
and relative pressure and kinetic energy difference between 
the levels; thus accurate measurement of these parameters 
is important for future research studies. 

High static stability is the only significant unique 
quality of the stratosphere available to explain the rela- 
tively slow development and large size of stratospheric 
systems compared to those of the troposphere. In de- 
veloping, these systems carry the seeds of their own 
destruction, for baroclinicity is dispersed and destroyed 
as the cold pool is transported from the darkness of the 
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Arctic into the sunlit latitudes while warm air floods the 
polar regions. Thereafter the long slow process of forma- 
tion of a new cold pool in the Arctic must begin again and 
continue until interrupted by the development of another 
perturbation or by radiational heating from the returning 
sun. The upper limit of the layer in which this process 
takes place has been the subject of much speculation. 
Strategically planned rocket firings near the poles in 
winter will certainly help to answer this problem. 

Much about the nature of the circulation in the upper 
stratosphere and lower mesosphere has been revealed by 
these and other rocket soundings. However, if anything, 
the nature of the circulation in the upper mesosphere 
has become more of a mystery. Even though diurnal 
and semi-diurnal effects and errors in measurement are 
large near the mesopause, the data presented here indicate 
that the circulation of nearby layers is really quite active, 
particularly so in winter through a deep layer at the 
boundary of the polar darkness, 
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1. INTRODUCTION 


In the eastern half of the United States the first half of 
January 1961 was relatively mild under high index circula- 
tion, continuing a warming trend which had started late in 
December [1]. During the second half of January, on the 
other hand, extremely cold weather prevailed in strong 
northerly flow between an intense trough along the east 
coast and a blocking High over western Canada; the cold 
spell lasted well into February. The highlight of the 
change in regime was a severe storm on January 19 and 20, 
1961, which caused heavy snow in major cities along the 
east coast from Washington, D.C. to Boston, Mass. Se- 
vere blizzard conditions were reported in southern New 
England where northeasterly gales caused extensive flood- 
ing of coastal lowlands. 

At no time during the month was the circulation favor- 
able for extensive precipitation over most of the interior of 
the contiguous United States. This was the driest January 
in 90 vears or more at several stations listed in table 1. 


2. TRANSITION WITHIN THE MONTH 


Many features of the interesting evolution of the mean 
circulation can be seen in the charts of 5-day mean 700-mb. 
height and weekly surface temperature anomalies, figures 
1 through 4. During the first week the mean circulation 
January 3-7, fig. 1A) consisted of a deep trough from 
Alaska southward, a rather strong ridge over the north- 
western States, and a trough along the east coast of North 
America from Labrador to Florida. A blocking ridge in the 
Greenland area, with largest height anomaly in Davis 


Tan_e 1,.—Stations reporting long-standing records for least January 
precipitation broken in 1961 


Date 


records 
began 


Station 


Grand Junction, Colo__. } 1892 
Detroit, Mich........-..- 
Escanaba, Mich-__. on 

Duluth, Minn_-___.._--- 

Rochester, Minn____.-_-- ; 

Kansas City, Mo.*___.. 

Akron, Ohio______. e 
Sandusky, Ohio________- 

Toledo, Ohio____. 

Chattanooga, Tenn_- 





*Equaled record. 


Strait, was joined to the subtropical Atlantic ridge, and a 
strong cyclonic center was located over the North Sea. 
The zonal index (not shown) was diminishing slowly from 
the early January peak which terminated the index cycle 
of December [1]. However the warming trend noted at 
the end of December over the contiguous United States 
continued, and temperatures over half or more of the 
country averaged above normal (fig. 1B). 

By the middle of the second week (January 10-14, fig. 
2A), the low center south of Alaska had intensified, while 
a blocking surge (best seen in the height anomaly pattern) 
had moved southwestward from Davis Strait to the Great 
Lakes. A weak trough over the Central Plains moved 
slowly eastward, and the trough in the western Atlantic 
grew stronger in the north. While zonal index values con- 
tinued to diminish slowly, they remained above normal. 

This was the warmest week of the month over most of 
the Nation (fig. 2B). Temperatures averaged above nor- 
mal over almost the entire country, averaging 10° to 25° 
above in the north central interior. 

One major storm during the week was responsible for 
moderate to heavy snow from West Virginia to southern 
New England. Part of this storm intensified off the mid- 
dle Atlantic coast, causing hurricane-force winds and high 
seas off the New England coast. 

Rapid amplification occurred from the central Pacific 
to eastern North America during the third week. Mid- 
way through the week, the Pacific cyclone reached its 
greatest intensity (over 1000 ft. below normal, fig. 3A). 
In the ridge downstream, heights averaged some 800 feet 
above normal in western Canada. With the Atlantic 
coastal trough also strong, the flow over central and 
eastern Canada was meridional, and cold air flowed into 
the eastern United States. 

These conditions accompanied the severe storm of the 
19th, which deepened more than 1 mb. per hour for 24 
hours over the relatively warm water along the Atlantic 
coast. This development in turn sustained and intensi- 
fied the meridional flow, so that the warm temperature 


regime of previous weeks was effectively ended in the 


East. 
22 were established in the wake of the storm at 


New record low temperatures for January 21 and 


New 


* Articles describing the weather of February, March, and April 1961 will appear in the 
May, June, and July 1961 issues, respectively. 
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Figure 1—(A) 5-day mean 700-mb. heights (solid) and departures 
from normal (dotted), both in tens of feet, for January 3-7, 1961. 
This represents the mean circulation during the first week of 
January. Note the positive anomaly center over Davis Strait, 
indicating the presence of blocking. (B) Departure of average 
temperature from normal (°F.) for the week ending January 8. 
Continuation of warming from the previous week was the general 
rule, except in the South (from [6)]). 


Haven, Conn., Wilmington, Del., Philadelphia, Pa., and 
for the 22d at Hartford, Conn., Chattanooga, Tenn. and 
Davtona Beach, Fla. Snowfall ranged from 2 to 16 inches 
in Maryland, 10 to 30 inches in New Jersey, and 10 to 20 
inches over southern New England. 

Intense cold dominated the weather east of the Con- 
tinental Divide during the final week. The long-wave 
pattern over North America (fig. 4A) remained strongly 
amplified, and northerly winds continued to transport 
frigid air into the United States. Sub-zero minima were 
reported southward as far as Denver, Colo., St. Louis, 
Mo., Nashville, Tenn., and Richmond, Va. In the South, 
only Brownsville, Tex., and southern Florida escaped 


freezing temperatures. Most of the snow cover over the 


Figure 2.—(A) 5-day mean 700-mb. heights (solid) and depar- 
tures from normal (dotted), both in tens of feet, for January 
10-14, 1961. By this time a positive anomaly center had 
retrograded from Davis Strait to the Great Lakes region, and 
the negative anomaly center south of Alaska had _ intensified 

(B) Departure from normal of average 

This 


from the previous week. 
temperature (°F.) for the week ending January 15, 1961. 
was the warmest week in January (from [6]). 


Northeast remained from the previous week, and addi- 

tional snow accumulated during a final storm which 

intensified off the Atlantic coast on the 27th. 
Comparison of half-month mean charts at 700 mb. in 


figure 5 reveals the remarkable change in circulation which 


occurred from the central Pacifie to the east coast of 
North America during the month. Contrast between the 
relatively zonal pattern of the first half (part A), and 
the strongly amplified pattern of the second half (part B) 
was greatest in the Yukon. 
(part () were as much as 800 feet 


Positive height changes there 
, While changes in the 
opposite sense were about half this magnitude in the 
troughs south of the Aleutians and over the Great Lakes. 


Changes of 700-mb. height anomaly were in the direction 
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Figure 3.—(A) 5-day mean 700-mb. heights (solid) and departures 
from normal (dotted), both in tens of feet, for January 17-21, 
1961. By this time the mean wave was strongly amplified from 
the eastern Pacific to the western Atlantic. An intense storm 
developed off the middle Atlantic coast this week and tempera- 
tures tumbled in the East as shown in (B) Departure from 
normal of average temperature (°F.) for the week ending January 
22, 1961. Most of the cooling occurred after mid-week (from [6]). 


favoring lower temperatures east of the Rockies. The 
close relationship of positive height anomaly in the Yukon 
and negative temperature anomaly in eastern United 
States is well known |2,3]. The change of half-month 
temperature anomaly (fig. 5D) shows that cooling occurred 
over the entire area east of the Continental Divide and 
Was most intense over Illinois and Indiana. 

During the latter half of January a 5-day mean High 
traversed an unusual path from the Great Basin north- 
ward along the western slopes of the Rocky Mountains 
through Alaska. The track of the High center has been 
plotted on the 15-day mean 700-mb. chart for January 
16 30 (fig. 5B), during which period the mean ridge over 
the Rockies was exceptionally strong. It is noteworthy 


Figure 4.—(A) 5-day mean 700-mb. heights (solid) and departures 
from normal (dotted), both in tens of feet, for January 24-28, 
1961. Large amplitude of the circulation continued from the 
previous week over North America. Cold air flooded the United 
States eastward from the Continental Divide. (B) Departure 
from normal of average temperature (°F.) for the week ending 
January 29, 1961. This was the coldest week of the month 
(from [6]). 


that the center followed closely the sharp ridge, and the 
High could be traced clearly on 5-day mean maps (pre- 
pared three times a week) from January 14 to February 4. 


3. MONTHLY CIRCULATION AND TEMPERATURE 
PATTERNS 


The most intense anomaly of the 30-day mean circula- 
tion at 700 mb. was centered in the extensive negative 
area in the eastern Pacific (fig. 6), associated with south- 
eastward displacement from normal of the Aleutian Low. 
In the center, about 700 miles south of the Alaskan 
Peninsula, 700-mb. heights averaged 470 feet below nor- 
mal and sea level pressures (fig. 7) 20 mb. below normal, 
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Figure 5.—(A) Mean 700-mb. contours (tens of feet) for January 1-15, 1961 and (B) January 16-30, 1961. Note amplification from 
(A) to (B) and path of a 5-day mean High over western North America from January 14 to February 4 in (B). Middle day of period 
is plotted above position and central intensity (in tens of feet) below. (C) Changes of 700-mb. height (tens of feet) (January 16-30 
minus 1-15). Rises over the Yukon contributed to falling temperature anomalies in the United States shown in (D), the tempera- 
ture anomaly change (°F.) for same periods as (C). 


Progressively lesser anomalies appeared downstream, in A secondary maximum of mean wind speed from the 
the ridge favored by climatology over the northern Arctic Basin southward through western Canada marked 
Rockies, in the trough along the east coast from the approximately the principal path of daily Highs (fig. 9A). 
Maritime Provinces southward, and in the Azores High. These Highs were instrumental in transporting cold air 
Blocking, as it is represented by positive height anom- into the eastern States in the latter half of the month. 
alies centered north of negative anomaly centers or During the warm portion, until around the 19th, all the 
channels, occurred in the Davis Strait area and over the  Aretic Highs moved eastward between Hudson Bay and 
East Siberian and Beaufort Seas. The Davis Strait region the Great Lakes. Subsequently the Highs plunged south- 
has been the seat of positive anomaly centers in most ward across the Great Plains and then followed the mean 
Januarys since 1954. Blocking is implicit in the spread- axis of maximum wind speed (fig. 8) eastward into the 
ing of the mean contours over Europe, though the anomaly Atlantic. In spite of the abnormality of the upper-level 
pattern there indicates that the condition was neither well- — circulation, it is noteworthy that both of these anticyelone 
organized nor strong. However, the axis of maximum paths are indicated as principal tracks for January by 
average wind speed at the 700-mb. level was diverted Klein [4}. 
southward from its usual path into central Europe (fig. 8). Although there was considerable variability within the 
Southward displacement also occurred in the eastern month, temperature anomalies over the contiguous United 
Pacific, eastern United States, and western Atlantic States (fig. 10) did not change greatly from Decem)er. 
Ocean. The usual test for persistence (Namias [5]) showed 85 
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Ficure 6.—Mean 700-mb. contours (solid) and departures from normal (dotted), both in tens of feet, for January 1961. Greatest anomaly 


was located in the eastern Pacific. 


percent zero or one class change (out of five classes), which 
is high compared to the expected 71 percent. During 
January temperature persistence was highest over the 
West and the South, where changes in the circulation 
(fig. 5C) were not extensive. One new record for Janu- 
ary cold was established, at Hartford, Conn., and several 
hew marks were set for warmth in southern California. 
Stations listing record high averages were Blue Canyon, 
Burbank, Los Angeles (city office), and San Diego. These 
Stations were in the region of warm downslope winds, 
locally called Santa Ana, which blow when the Great 


Ridges and troughs were progressively weaker going downstream. 


Basin High is strong at sea level. The mean sea level 
pressure distribution of January (fig. 7), with pressures 5 
mb. higher than normal over the Great Basin, was ideal 
for the generation of these warm winds. In other western 
regions temperatures were generally warmer than normal 
under the influence of the persistent mean ridge. 


4. PRECIPITATION 


It has been mentioned that the interior of the contiguous 
United States was extremely dry in January. 


Figure 11, 
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Figure 7. 
normal (dotted), (both in millibars) for January 1961. 
mediate isobars are dashed. 
normal south of Alaska, 
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Figure 8.— Mean isotachs (in meters per second) of 700-mb. wind 
speed during January 1961. Solid arrows indicate observed axis 
of maximum wind and dashed arrows the normal. 
Secondary maximum approximated principal path of Arctic High 
from the Polar Basin across western Canada into central United 
States (fig. 9A). 


speed, 








Mean sea level isobars (solid) and departures from 
Inter- 
Pressures were 20 mb. lower than 











B. CYCLONES TRACKS 2.3 2 OUTS 
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FIGURE 9. 
cyclone passages (within equal area boxes of 66,000 nautical 
mi.”) during January 1961. Aretic Highs were diverted toward 
James Bay during first 20 days of month, but continued south- 
ward through Great Plains thereafter. Areas of zero frequency 
are shaded, 


(A) Frequency of daily anticyleone passages and (B) 


showing the percent of normal monthly precipitation, 
illustrates the large extent of the area reporting less than 
normal amounts. In addition to the broken records of 
long standing listed in table 1, there were numerous reports 
of “driest January” for shorter periods. Apparently the 
strong ridge over western North America precluded ex- 
tensive upward motion and flow of appreciable moisture 
from the Gulf of Mexico or Pacific Ocean. Certainly no 
daily cyclones passed through the ridge at sea level (see 
fig. 9B). Several Lows of the Alberta type migrated 
eastward into the Great Lakes region, but Alberta Lows 
are typically dry. 

Another group of sea level cyclones originated in Kansas 
or the Texas Panhandle. 


These too were dry initially, 
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Ficure 10.—Departure of average temperature (°F.) from normal 
for January 1961 (from [6]). Hatching indicates areas of above 
normal temperatures. 


in the lee of the Divide, but produced precipitation farther 
eastward when moisture entered the systems from the Gulf 
of Mexico. This group of Lows triggered the intense 
developments off the Atlantic coast described in section 2. 
Precipitation totals in the region affected by these storms 
averaged Jess than the January normal amounts except 
for a small area in the central Appalachians and a narrow 
strip along the coast from Cape Hatteras to Long Island. 
Snowfall was greater than normal over a larger area. 
Harrisburg, Pa., reported 34 inches, a new January record. 
At Akron, Ohio, snowfall was twice the normal, but 
its water content was very low and total monthly pre- 
cipation was the lowest on record. 

In central Texas precipitation amounts were as much 
as 400 percent of normal. Abilene and Waco reported the 
wettest January on record with 3.99 inches and 5.83 inches, 
respectively. There the precipitation was associated with 
over-running Gulf air during the passage of upper-level 
Lows or troughs across the Southwest, none of which was 
reflected appreciably in the daily sea Jevel pressure 
patterns over Texas. 


Fiavure 11.—Percentage of normal precipitation for January 1961, 


6. 


The month was extremely dry over most of the interior of the 


United States (from [6]). 
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Weather Bureau Technical News 


SPECIAL SLIDE RULE FOR COMPUTING EVAPOTRANSPIRATION 


In response to the need for a rapid and accurate method of 
making daily computations of evapotranspiration from forage 
crops, a special slide rule has been developed by L. T. Pierce, 
Weather Bureau State Climatologist in Columbus, Ohio, and G. H. 
Snyder, Student Trainee in 1959 and 1960 at the Weather Bureau 
Office in Columbus. Use of this special rule, which incorporates 
a moderately complex empirical formula, has cut the time consumed 
in making numerous daily calculations by at least 50 percent of 
that formerly required when using a conventional slide rule or 
desk calculator. At the same time it has reduced the chance of 
error and eye strain from constant reference to tables. 

This particular slide rule was designed for use in computing 
“actual”? evapotranspiration from established meadow crops in 
Ohio and adjacent areas through use only of the commonly available 
daily temperature and rainfall data. It is intended for current 
use as a device for keeping track of changes in soil moisture, and 
also for reconstructing the probable moisture conditions prevailing 
An adaptation of the same type of empirical formula 


in past years. 


is now being developed for use with other common crops such as 
corn and wheat; and computations for these crops can be made on 
the same slide rule by substituting appropriate scales. Other 
workers in meteorology and related fields might well find this kind 
of tool very useful for making routine computations involving 
empirical formulae. In several respects, the special slide rule is 
preferable to the commonly-used nomogram, especially when 
large numbers of computations are required. 

A full description of this rule, with particular emphasis on the 
procedures employed in its development, is given in a paper by 
L. T. Pierce and G. H. Snyder, ‘‘Development and Use of a Special 
Slide Rule for Computing Meadow Evapotranspiration,’’ which is 
available as a Weather Bureau “printed manuscript’’ from the 
Publications Section, U.S. Weather Bureau, Washington 25, D.C. 
This paper, which contains full-size scales for use in constructing 
a duplicate of the original rule, is a ‘“‘how-to-do-it’’ type of pres- 
entation that makes it relatively easy to adapt the procedures 
to other uses. 
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Suggestions for Authors 


Articles are accepted for the Monthly Weather Review 
with the understanding that they have not been published 
or accepted for publication elsewhere. 

Two copies of the manuscript should be submitted. All 
copy, including footnotes, references, tables, and captions 
for figures, should be double spaced with margins of at 
least 1 inch on sides, top, and bottom. Some inked 
corrections are acceptable, but pages with major changes 
should be retyped. The style of capitalization, abbrevia- 
tion, ete., used in the Review is governed by the rules set 
down in the Government Printing Office Style Manual. 

An abstract should be supplied with all but the briefest 
of papers. It should cover, in as concise a manner as 
possible, the main purpose and subject matter of the paper 
and indicate briefly the principal conclusions. One 
paragraph is usually sufficient. 

Tabies should be typed, each on a separate page, with 
atitle provided. They should be numbered consecutively 
in arabic numerals. 

In equations conventional symbols in accordance with 
the American Standards Association Letter Symbols for 
Meteorology should be used. If equations are written 
into the manuscript in longhand, dubious-looking symbols 
should be identified with a penciled note. 

References should be listed on a separate sheet numbered 
in the order in which they occur in the text; or, if there 
are more than 10, in alphabetical order according to 
author. The listing should include author, title, and 
source (if the source is a magazine, the volume, number, 
month, vear, and complete page numbers should be listed; 
if the source is a book, the publisher, place and date of 
publication, and the page numbers of the reference). 
If reference is made to a self-contained publication, the 
author, title, publisher, place of publication, and date 


should be given. Within the text, references should be 
indicated by arabic numbers in brackets to correspond to 
the numbered list. 

Footnotes should be numbered consecutively in arabic 
numerals and indicated in the text by superscripts. Each 
should be typed at the bottom of the page on which the 
footnote reference occurs. 

Illustrations.—A list of captions for the illustrations 
should be typed (double spaced) on a separate sheet. 
Every illustration should have an explanatory caption. 
Number each illustration in the margin or on the back 
outside the image area. To fit into the Review page, 
illustrations must take a reduction to 3% by 9 inches 
(column size) or 7% by 9 inches (page size). Map bases 
should show only political and continental boundaries 
and latitude and longitude lines, unless data are to be 
plotted, when station circles will also be needed. Usually 
the less unnecessary detail in the background the better 
will be the result from the standpoint of clear reproduction. 

Line drawings and graphs should also be uncluttered 
with fine background grids unless the graph demands very 
close reading. Lines should be even and black, lettering 
legible, symbols distinctive, and both lettering and 
symbols large enough to maintain their legibility under 
the necessary reduction. If reproductions of line drawings 
are submitted because of the large size of the original, 
they must be extremely clear. Mimeographed copies or 
“photocopy” types of reproduction are not satisfactory. 

Photographs should be sharp and clear, with a glossy 
surface. Bear in mind that marks from paper clips or 
writing across the back will show up in the reproduction. 
Drawings and photographs should be protected with 


cardboard and mailed flat. 
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Publication of Weather Notes 


Many years ago the Monthly Weather Review published detailed eye- 
witness accounts of exceptional storms. These accounts both enrich the 
meteorologist’s knowledge of storms and provide him with particular details 
that cannot be found elsewhere. Because such information bears directly 
upon questions the meteorologist must attempt to answer about weather 
phenomena (for example, the identification of storms as tornadoes), and 
because the information has potential value in both the research and service 
programs of the Weather Bureau, publication of eye-witness accounts and 
brief analyses of exceptional storms and other meteorological phenomena 
was resumed in the April 1955 issue. They appear from time to time under 
the heading “Weather Notes’’. 

Contributions to these ‘“‘Notes’’ are invited from readers of the Review. 
There is no limitation placed on length of description but it is expected that 
most will be short accounts. Any weather peculiarities, whether storms or 
other phenomena, are acceptable subject matter. Material should be addressed 
to Editor, Monthly Weather Review, U.S. Weather Bureau, Washington 25, 
D.C. 
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